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\q Abstract 

We compute the expectation value of the circular Wilson loop in Af = 2 
supersymmetric Yang-Mills theory with Nf = 2N hypermultiplets. Our 
O results indicate that the string tension in the dual string theory scales as 

the logarithm of the 't Hooft coupling. 
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1 Introduction and Summary 



Exact results in quantum field theory usually rely on powerful symmetry 
principles, such as supersymmetry. A beautiful example of non-perturbative 
use of supersymmetry is Pestun's exact calculation of circular Wilson loops 
in a wide class of M = 2 supersymmetric Yang-Mills theories [TJ . To be more 
precise, Pestun reduced the problem to a finite-dimensional matrix integral, 
still to be evaluated. In this paper we address this problem in the 't Hooft 
large- N limit. 

The large- N limit of Pestun's matrix model has been discussed in [2J, 
with rather unexpected conclusions on its strong-coupling behavior. The 
conclusions of [2J for the most part rely on simple scaling arguments. Our 
goal will be to develop systematic strong and weak coupling expansions of 
Pestun's matrix model and its observables. We shall see that the strong- 
coupling behavior of the matrix integral is indeed very unusual. The large- 
N master field has an infinite support at infinite coupling, but such that 
the majority of matrix eigenvalues remain finite in accord with the scaling 
arguments. The limiting master field describes a certain class of observables, 
but not all observables. In particular, it does not describe Wilson loops, 
which are determined by the largest eigenvalue as in ordinary matrix models. 

We shall concentrate on the M = 2 superconformal Yang-Mills theory 
(SCYM), an SU(N) gauge theory with Nf = 2N hypermultiplets in the 
fundamental representation. This theory has zero beta function and therefore 
is superconformal at any value of the Yang-Mills coupling g. One may expect 
that at iV — > oo and when the 't Hooft coupling A = g 2 N is large the 
SCYM theory is described by a weakly-coupled dual string theory. There is 
a number of proposals for the dual string/ supergravity background [3J IU E] ■ 
Neither of these backgrounds is of the simple form AdS^ x X 5 . In [31 [5], 
the AdS$ part of the type-IIA/M-theory geometry is warped with respect 
to the coordinates of X 5 . In [I], the string dual is non-critical, so X 5 is 
less than five-dimensional, or even partly non-geometric. In either case the 
dual string theory is never completely classical, since the backgrounds of 
[31 [5] contain a curvature singularity, and the non-critical string of pE] is 
intrinsically quantum. 

We believe that Wilson loop calculations (along with possible integrability 
of the SCYM theory [H]) may shed more light on the AdS/CFT duality in 
the M = 2 setting. For one thing, Wilson loops are very sensitive probes 
of the dual string dynamics, as they couple directly to the string worldsheet 
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[7J. One can also do exact field-theory calculation at any coupling with the 
help of Pestun's results [lj. In the more familiar M = 4 case, the exact 
calculation of the circular Wilson loop [H |9], [I] immediately confirms the 
AdS / CFT relationship between the string tension and the 't Hooft coupling 

T^ 4 =^. (1.1) 

Our goal will be to derive a similar relationship for M = 2 SCYM. 

The circular Wilson loops in M = 4 SYM were first calculated by re- 
summing planar diagrams (8J [9] . The planar perturbation theory for Wilson 
loops in Af = 2 SCYM has been studied in [TT] . where quite interesting reg- 
ularities have been observed. To make connection to this work we will also 
study the weak-coupling expansion of Pestun's matrix model. 

The paper is organized as follows: in the next subsection we briefly sum- 
marize our strong-coupling results and discuss their possible implication for 
the string dual of Af = 2 SCYM. In sec. [2] we review Pestun's matrix model 
for Af = 2 SCYM on S 4 and derive the saddle-point equations for the eigen- 
value distribution in its large- N limit. In sec. [3] we study the weak and in 
sec. [4] the strong coupling expansions of the saddle-point equations, and in 
sec. [5] we compute the one-instanton correction. We conclude with general 
discussion in sec. M The technical details of our calculations are collected in 
the appendices. 

1.1 Summary of results 

The vacuum expectation value of the Wilson loop is computed holographi- 
cally by summing over random surfaces in the bulk which end on the given 
contour on the boundary [TJ. Assuming that the string tension is large, the 
string path integral is saturated by the surface of the minimal area. For the 
circular loop on the boundary of AdS$ the regularized minimal area (which 
is negative for any contour) is equal to — In p2]. The expected form of the 
circular Wilson loop expectation value in the semiclassical regime then has 
the form 

W(C cilclc ) = K ■ T- 3 / 2 e 2 ^ T (T oo), (1.2) 

J By the string tension we will always mean the dimensionless ratio of the AdS radius 
squared to a': T = R 2 /2na' . 
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where T is the string tension. The factor of T~ 3 / 2 comes from the gauge 
fixing in the string path integral, as explained in {9j , and the constant factor 
K is determined by the quantum fluctuations of the string worldsheet. 

The semiclassical calculation of the circular Wilson loop was sensitive only 
to the universal AdS 5 factor in the geometry, that has to be there because of 
the conformal symmetry. The structure of the internal space {X 5 ) was not 
very important, and could only enter through the prefactor in the formula 



(1.2). In this respect the result (1.2) looks completely universal and should 
apply to any theory with the AdS dual. However, the derivation also assumes 
that the semiclassical approximation is accurate in some range of parameters, 
which might or might not be the case. It might happen that the string dual 
of = 2 SCYM is always in the quantum regime, either because of the 
curvature singularities present in the background geometry j3[ [5] , or because 
of strong quantum fluctuations that cancel the central charge in case the 
string dual is non-critical [I] . Whether or not the quantum sector in the string 
sigma-model can be separated from the geometric AdS$ factor in the Wilson 
loop calculation, and whether the curvature of AdS$ ever becomes small or 



not is unclear to us. In some sense ( 1.2 ) can be regarded as a parameterization 



of the Wilson loop vev in terms of the effective string tension. We will find 
that the circular Wilson loop vev in M = 2 SCYM is indeed consistent with 



the parameterization (1.2) upon a simple identification 



T=ll„A. (1.3) 



This is certainly very different from the standard M = 4 relation (1.1). 

The logarithmic behavior of the effective string tension follows from the 
power-like growth of the Wilson loop expectation value at strong coupling: 

A 3 

W(C circlc ) = const — . (1.4) 

(In A) 1 

We will try to be accurate with normalization, although keeping constants 
on top of the logarithms is always difficult. Our estimate for the coefficient 



K in (jl_2|) is 

if ~ 3.13 • 1(T 5 . (1.5) 



The corresponding constant in (1.4) is 9.47 • 10 . These are analytic but 



approximate predictions, which we expect to have a few percent accuracy. 
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2 Partition function and saddle-point equa- 
tions 



The field content of the Af = 2 SCYM theory consists of the SU(N) gauge 
field two adjoint scalars $/ from the vector multiplet (I = 1,2), 2 x 
2N fundamental scalars Qf, Q^ A from 2N hypermultiplets (A = 1,2; / = 
1, . . . ,2N), and various fermions that make the spectrum supersymmetric. 
The Wilson loop is defined as 



W{C) 



1 

N 



tr P exp 



ds (iA^x)^ + n I ^ I (x)\x\) 



c 



(2.1) 



where n is a unit two-dimensional vector, which can vary along the contour, 
but in Pestun's calculation has fixed constant components. 



2.1 Partition function on S 4 

Pestun computed the partition function of M = 2 SCYM on S A using localiza- 
tion p. Compactification on the sphere provides a useful IR regularization, 
but otherwise is not important, since we are dealing with a conformal theory 
and the sphere is conformally equivalent to M 4 . Let us briefly review the 
main steps of Pestun's calculation. 

One can define a "vacuum state" of M = 4 SCYM "on the Coulomb 
branch" by Higgsing the theory with an expectation value of the scalar field 
$/ along the n direction. The expectation value can be brought to the 
diagonal form, with respect to the SU(N) indices, by a gauge transformation: 
($/) = 7ifdiag(ai, . . . , a^). The eigenvalues a, should satisfy the SU(N) 
constraint: 

5^ = 0. (2.2) 

i=l 

On the sphere one actually has to average over all vacua, but this can be 
postponed till the last moment. A useful strategy is to first compute an 
effective action for aj by integrating out all other fields. Because of the su- 
persymmetry this effective action is one-loop exact. A proof and the explicit 
calculation of the one-loop factors can be found in pp. 
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The path integral of M = 2 SCYM on S 4 thus reduces to a matrix integral 
over the zero mode of the adjoint scalaiF] 



Z = d N 1 a Y[ (a* ~ «j) 2 e ~^ ' 1 Z 1Aoop (a) \Z inst (a; g 2 ) \ 2 (2.3) 



The first factor in the integrand (the Vandermonde determinant) is the 
Faddeev-Popov determinant of the diagonal gauge. The exponent in the 
second factor is the classical action - the area of the four-sphere times the 
conformal coupling of the scalar to the constant curvature of S 4 . The last 
two terms are the one-loop determinant from integrating out field fluctua- 
tions and the instanton contribution. 

The instanton partition function Zi nst (a; g 2 ) is a known |13j . albeit fairly 
complicated function of the eigenvalues and the Yang-Mills coupling. It 
is usually assumed that instantons are not important in the large- N limit, 
because the instanton action is linear in N in the large- iV limit: 




In practice the situation is more complicated, because the instanton weight 
contains moduli integration and the number of instanton moduli grows with 
N . Whether instantons are suppressed or not is thus a dynamical ques- 
tion |14j . There are known examples where the volume of the moduli space 
overcomes the suppression by the instanton action and leads to a large-iV 
phase transition into a non-perturbative phase. Later we will compute the 
one-instanton weight to check if there is an instanton-driven phase transi- 
tion in M = 2 SCYM. We will find that the moduli integration enhances 
the instanton weight by a factor of \/N, which is insufficient to overcome 
the exponential suppression by the instanton action. For now on we just set 
Z inst (a;g 2 ) = 1. 

The scalar vev acts on the adjoint fields from the vector multiplet through 
the commutator: [(<&) ,v]ij = (d{ — dj)vij and therefore the vector multiplet 
fields get masses mf,- ~ (cij — cij) 2 - On the fundamental fields in the hy- 
permultiplets the scalar vev acts by multiplication: (($) h)i = Oj/ij and the 
hypermultiplets get masses m 2 ~ a 2 . The one- loop effective action is thus 

2 We choose to normalize the Yang-Mills kinetic term as — tr F^ l/ /2g 2 , which differs by 
a factor of 2 from the conventions in pQ . 
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a combination of terms that depend on 



) 2 and a\. The explicit ex- 



pression was computed by Pestun, and can be expressed through a single 
function H(z), which is related to the Barnes G-function: 



H(x) 



-(1+7)^ 



G(l +ix)G(l -ix), 



and admits an infinite-product representation: 

oo 

//(*•) = n ( 1 + 



n=l 



X 



(2.4) 



The one-loop factor in the partition function is given by [T] 



loop 



li< 3 



(2.5) 



where the numerator is the contribution of the vector multiplet and the 
denominator is the contribution of 2N hypermultiplets. 



The expression (2.5) has an interesting symmetry. If we multiply H(z) 
by a Gaussian: 



H(z) H(z)e 



Cz 2 



(2.6) 



with an arbitrary constant C, the partition function (2.5) will not change 



because of the trace condition (2.2). This property is closely related to 
the UV finiteness of M = 2 SCYM. Each individual one- loop determinant 



in the background field ($) diverges in the UV. In (2.4) this divergence is 



regularized in an arbitrary way, but in a more general setting, for instance 
in the M = 2 theory with Nf ^ 2N, the coefficient in front of a 2 in the 
effective action will be logarithmically divergent and will require adding a 
counterterm that shifts 1/A and results in the one-loop beta-function. In 
M = 2 SCYM, any term quadratic in a cancels due to the shift symmetry 



(2.6), which in particular means that the beta function is zero. It also means 
that there are no finite corrections to the coefficient in front of a 2 . The 
absence of finite renormalization has interesting consequences for the Wilson 
loop vev. The first renormalized vertex in the effective action is thus quartic. 
The lowest-order diagram with a quartic vertex contains three propagators 
(fig. [I]) and thus contributes to the Wilson loop vev at 0(A 3 ). The lower 
orders are described by the Gaussian matrix model and thus are the same as 
in the M = 4 super- Yang-Mills, in agreement with the analysis of the SCYM 
perturbation theory 
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Figure 1: The lowest-order non-Gaussian correction to the Wilson loop in the 
matrix model appears at order 0(A 3 ). 



The circular Wilson loop in SCYM on S* 4 is a chiral observable and there- 
fore is directly related to Wilson loops (matrix exponentials) in the matrix 



model (2.3): 



W(C cir , 



irclc I 



(2.7) 



2.2 Saddle-point equation 

In the large-iV limit the saddle-point approximation becomes exact for the 



integral (2.3). The effective action for the eigenvalues is 



S(a) = J2( S Y^ + 2\nH(a i )] - ^ ^ (in (a* - a 3 ) 2 + 21n#(a* - a,)) . 

i ' i<j 

(2.8) 



Minimizing the action we get the saddle-point equations: 

r 2 



57T 

T 



o< - Kfa) - — ( — K(ai- aj) ) = 0. 



(2.9) 



The function K(x) that appears here is defined as 

H'(x) 



K{x) 



H(x) 



(2.10) 



and is given by 



oo , 
n=l ^ 



Tl 

' x (ip (1 + ix) + ip (1 - ix) - 2^(1)) 



(2.11) 
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where ip(x) is the logarithmic derivative of the Gamma-function: ijj(x) 



T'(x)/T(x). The symmetry (2.6) translates into invariance of the saddle- 
point equations under the shifts of K(x) by a linear function. We have used 
this symmetry to normalize K(x) such that its Taylor expansion starts with 
0(x 3 ). 

As usual, the saddle-point equation can be interpreted as an equilibrium 
condition for N pairwise interacting particles in the common external poten- 
tial. Since 

K(x) « 2C(3)x 3 {x->0), (2.12) 

the K(x) terms in the force are negligible at short distances, and the equi- 
librium distribution is determined by the balance of the attractive harmonic 
potential and the 1/x pairwise repulsion. At large distances, on the contrary, 
K(x) is large: 

K{x)^2x\nx (x^+oo), (2.13) 

so the total one-body potential becomes repulsive, while the two-body inter- 
action becomes attractive, exactly opposite to the short-distance behavior. 
The balance with this configuration of forces is still possible, but the sys- 
tem becomes potentially unstable towards the spread of the eigenvalues to 
infinity. The stability again follows from mutual cancellation between the 
one-body K{a,j) and two-body K(cii — aj) terms for each individual particle, 
such that at large distances the harmonic potential still plays the most im- 
portant role and confines the eigenvalue distribution to a finite interval. The 
cancelations have basically the same origin as the UV finiteness of theory. 
We can introduce the eigenvalue density, 

P( x ) = jfY, 5 (2.14) 

i 

which is defined on some interval (— and is unit normalized. The saddle 
point equations take the form of a singular integral equation: 

dyp(y)(—-K(x-y)) = 8 ^x-K(x), (2.15) 



x — y J A 



which together with the normalization condition determines the density and 
the endpoint \i. The Wilson loop expectation value is given by the Laplace 
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transform of the density: 



circle J 



dx p{x) e 



2-kx 



(2.16) 



We will develop systematic expansions of the saddle-point equation (2.15) 



at weak and at strong coupling. To this end, it will prove useful to rewrite 
the equation in a different form, suggested by the solution of the Hermitean 
one-matrix model [15] (which corresponds to setting K(x) to zero). In this 



case the integral operator in (2.15) can be inverted by application of 



/ 



dx 



'p 2 -x 2 z-x 



to both sides of the equation. If we apply this operator to (2.15), we get: 



p(x) 



8tt 

T 



l 



dy 



' p 2 — x 2 



it—- /-v»2 

tx £ J x — y y p — y 

-A* 



dzp(z) (K(y-z)-K{y)). 

(2.17) 

We thus reduce the problem to an integral equation of the Fredholm type, 
since the kernel is not singular any more. The endpoints of the eigenvalue 
distribution are determined by the normalization condition: 



4vr 2 /i 2 1 

A 7T 



dyy 



vV -y 2 



dzp{z){K{y-z)-K{y)). (2.18) 



3 Weak coupling 



In the weak coupling regime when A <C 1, it follows immediately from the 



saddle point equation (2.15) that the eigenvalues are distributed on an inter- 
val (— p, p) with /i<C 1. In this regime, we can therefore express the function 
K(x) by the Taylor series 



K(x) = -2^(-l)"C(2n + l)x 



2n+l 



(3.1) 



n=l 
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where ((n) is the Riemann zeta function. At the lowest order of approxima- 
tion we truncate the series by keeping only the first term (2.12). 

The last term in the equation (2.17) then factorizes, the integrals over z 
and y can be done separately, and give an approximated expression for the 
density 



p(x) 



?7r 6C(3)m 2 
X + n 



where we defined the second moment 

1712 = J dz p(z)z Z 



(3.2) 



(3.3) 



The moment m 2 is self-consistently determined by the density (3.2), and we 
obtain the equation 



m 2 



87r 6C(3)m 2 ^\ vr/i 4 



(3.4) 



that permits us to express m 2 in terms of p and A. Therefore, at this or- 
der of approximation, the density and the normalization condition are given 
explicitly by 



p(x) 



8tt 24vrC(3)/i 4 



A 4A-3C(3)A/x 4 
4vrV 12vr 2 C(3)/i 6 



/ /i 2 — x 2 



A 



+ 



4A - 3C(3)A/i 4 



(3.5) 
(3.6) 



and the normalization condition can be used to express p as a function of A. 
It results in 

VX 3C(3)A 5 / 2 

V= — org-K +••• ( 3 - 7 ) 



2tt 



256tt e 



The weak coupling solution (3.5 ) is compared to the numerical data in fig. |2j 
Using the density (3.5) and the expression (3.7) it is possible to compute 
the expectation value of the circular Wilson loop, as described in (2.16). We 
obtain 

W(C circlc ) = 1 + A A " ( 1 



8 192 



9216 



3C(3) \ 3 , 
512^J A + 



(3.8) 



that is in agreement with the result of [TT], considering the large N limit. 
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0.000 



0.005 



0.010 



0.015 



Figure 2: The weak coupling density at A = 0.01 (solid line), compared to the 
numerical solution of the saddle-point equation with N = 1200, shown in circles. 

The perturbative scheme just outlined can be pushed to an arbitrary 
high order in A. In particular, considering an expansion for K(z) up to order 
0(z 2M+1 ), the Fredholm equation (2.17) gives an approximated expression 
for the density that depends on M moments m^i with % = 1, . . . , M. This 
density can be used to compute the moments, giving a system of M equations 



that is the higher order generalization of the (3.4) and that permits us to 



express the moments as functions of \i and A. In this way the density and the 
Wilson loop can be computed to an arbitrary order in A. Expanding K(z) 
up to order 0(z 2M+1 ), it is possible to compute the expectation value of the 
Wilson loop up to order 0(X 2+M ). In the appendix |A| this scheme is pushed 
to the seven-loop order for the Wilson loop vev. 



4 Strong coupling 

The strong-coupling behavior of M = 2 SCYM on S A was analyzed in [2]. 
Using simple scaling arguments the authors of [2] reached the conclusion that 
the eigenvalue density approaches a finite limiting shape at A — > oo. Let us 
denote this limiting shape by Poo{x). This behavior is in sharp contradistinc- 
tion to the Gaussian matrix model which describes Af = 4 super- Yang-Mills 
[8j [9] , where the density is a function of the scaling variable x / For the 
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M = 2 SCYM matrix model no consistent scaling is possible, as shown in 
[2], and therefore the density freezes out as A — > oo. On the other hand, the 



weight in the integral representation of the Wilson loop (2.16) exponentially 
grows with x, and thus the Wilson loop vev is determined by the largest 
possible eigenvalue, x — W ~ e 2nfJl [2J. In the Af = 4 context this leads 
to an exponential growth of the Wilson loop vev with \/A, and confirms the 



square-root dependence of the string tension on the 't Hooft coupling (1.1). 



The edge behavior of the density determines the prefactor in (1.2): the factor 
of T -3 / 2 arises as a consequence of the square-root singularity at the end- 
point of the eigenvalue distribution [2]. If in Af = 2 SCYM the endpoint 
approaches a finite limiting value p^, the Wilson loop vev also approaches a 
constant, which is really difficult to reconcile with the holographic duality. 
However, assuming that p,^ is a finite constant that does not depend on 



A we immediately run into contradiction. The integral equation (2.15) holds 
on the whole interval (— p, p) and in particular should be satisfied at x = p. 
Assuming that the density freezes out at A — > oo we can set A = oo in the 
equation and then at the endpoint we get 

dypoo{y) (k(hoo) - K{p OQ —y) + 



Poo - y, 

"Moo 

where we have used the normalization condition to move the first term inside 
the integral. The integral here converges without any regularization, because 
the density goes to zero at the enpoint as a square root of the distance: 
Poo(y) ~ v 1 Poo — y- But K(x) is a monotonically growing function, and 
hence the integrand is strictly positive, so the equality can never be satisfied. 

What can resolve this paradox? We must assume that the density ap- 
proaches a limiting shape at A = oo, otherwise we run into contradiction 
with the scaling arguments of [2] . The only way to reconcile the existence of 
the limiting density with the saddle-point equation is to admit that p^ = oo. 
This behavior is extremely unusual for matrix models. It is ultimately related 
to the intrinsic instability of balancing the one-body repulsion against the 
two-body attraction, as we discussed in sec. [2] We have performed extensive 



numerical checks of this behavior by solving the saddle-point equations (2.9) 
at large but finite N and 1/A = 0. In the numerics we have dealt with finitely 
many eigenvalues which of course do not extend to infinity, but the largest 
eigenvalue, as we have found, grows more or less linearly with N without any 
signs of saturation. 
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Figure 3: The density at infinite coupling (solid line), compared to the numerical 
solution of the saddle-point equation with N = 1200 and 1/A = 0, shown in circles. 



4.1 Infinite coupling 

The equation for the limiting shape of the eigenvalue distribution is 

+oo 

-j dypoo(y) (^T—y ~ K ( x ~ V)j = ~ K ( X )- 



(4.1) 



It can be solved by Fourier transform: 

1 



m sign u 



1 + 



2 sinh 2 



Poo(^) 



m sign oo 
2 sinh 2 ^ ' 



which gives: 



and 



Poo(^) 
Poo{x) = 



1 



cosho; 
1 

2 cosh if 



(4.2) 

(4.3) 
(4.4) 



In fig. [3] this solution is compared to the numerical data. 

The majority of eigenvalues are indeed concentrated at x = 0(1), as the 
scaling arguments suggest, but the density has exponential tails extending 
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all the way to infinity. The shape of the density is such that the scaling 
arguments apply to an observable (0(x)) if and only if 0{x) does not grow 
with x faster than e 7 ™/ 2 . The expectation value (O(x)) then approaches a 
constant limiting value at A — > oo. The Wilson loop does not belong to this 



class of observables, since the integral (2.16) with the limit density diverges 



Therefore, the strong-coupling behavior of the Wilson loop is still determined 
by the largest eigenvalue and is sensitive to the precise value of \i and to the 
endpoint behavior of the eigenvalue density at large but finite A. In order 
to compute the Wilson loop we need to analyze deviations from the strict 
infinite-coupling limit. Technically this is a much more difficult problem, 
and we start with simple qualitative estimates before proceeding to a more 
systematic analysis. 

4.2 Qualitative estimates 

If A is large but finite, the density cannot differ much from Poo(x). The most 
important difference is the finite extent of the eigenvalue distribution. The 
exponential tails of Poo{x) are trimmed at x — ±/x in a certain way, which we 
are going to analyze. As we have seen, there is an intrinsic instability in the 
balance of forces at strong coupling. Because of this instability the effective 
force acting on an individual eigenvalue goes to zero at infinity (when A is 
strictly infinite), and this allows the eigenvalues to spread all over the real 
line. If A is not infinite but very large, the linear term in the force can 
be neglected in a certain range of x, since it has a very small coefficient. 
But eventually it takes over and starts to press the eigenvalues towards the 
origin. As a result, the density becomes slightly bigger everywhere, and 
cannot extend further than x — zkfi. As a crude approximation we can 
neglect the linear force at x ~ 1 and neglect the two-body forces at x ~ /i. 



The eigenvalues density is then approximated by (4.4) at x ~ 1 and by the 
Wigner semi-circle at x ~ p: 

P{%) ~ y vV - x2 ■ ( 4 - 5 ) 

Matching the two at some intermediate scale we find that the endpoint fi 
should scale logarithmically with A: 

2 

/ u = -lnA + ... (A^oo). (4.6) 

7T 
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We can obtain a more accurate estimate by taking into account the nor- 
malization condition. The linear force eliminates all the eigenvalues that sit 
at x > fi, compressing them to smaller x. Let us consider the positive end- 
point x = +fi. The fraction of eigenvalues redistributed to smaller x roughly 
speaking is given by 

oo 

J dx Poo (x) ~ ^ e~^ 2 . 

If we approximate the eigenvalue distribution near the endpoint by the 
Wigner semi-circle, the excess number of eigenvalues, compared to the in- 
finite-range case, is obtained by integrating the Wigner distribution from 
some cutoff scale fi — zq (zq ~ 1) to fi: 

/g7T 10 7t Z^^^ 
dx — vV -x 2 ~ 3A ° v 7 ^ • 

Equating this with the number of eigenvalues that came from infinity, we get 
an equation 

Cy^e^^A, (4.7) 

which determines /i as a function of A. We cannot find the constant of pro- 
portionality in this equation from the simple arguments above (the constant 
depends on the cutoff zq which we have put in by hand). To honestly compute 
this constant we need more sophisticated methods. Later we will estimate 

C ~ 14.60, (4.8) 

which corresponds to taking Zq = 0.54 in the previous simple-minded argu- 
ment, from which we can get an idea how well the actual density near the 
endpoint is approximated by the Wigner semi-circle. 

In fig. [4] we compare the weak-coupling and strong-coupling predictions 
for the endpoint position with the numerical data (see also appendix [B| . 

We can estimate the Wilson loop vev in a similar wajj^j Here we may 
expect better accuracy, because the exponential weight of the Wilson loop is 
peaked near the positive endpoint of the eigenvalue distribution, and this is 

3 We would like to thank the anonymous referee of JHEP for suggesting this calculation 
to us. 
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Figure 4: The ratio /i/vA as a function of A: dots represent numerical results; 
black curve is the seven- loop weak-coupling approximation (A.7); purple curve is 
the strong-coupling prediction (4.7), (4.8). 



precisely where it makes sense to use the simple semi-circle approximation 



(4.5): 



W(C t 



circle J 



87T 

T 



(4.9) 



Doing the integral we find that the Wilson loop vev scales with A and \x as 

(4.10) 



W(C circlc ) = i?Y e2 ^ 



with R = 2. Later we will get a better estimate of the proportionality 
constant which, as a matter of fact, does not differ much from the simple- 
minded approximate calculation above: 



R = 2.18. 



(4.11) 



Expressing jj, in terms of A with the help of (4.7), we find from (4.10) the 



results (1.4), (1.2) quoted in the introduction. For the constant K in (1.2) 
we get: 



K 



3/2 



R 



(4.12) 
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and using (4.8), (4.11) we obtain the numerical estimate (1.5). 

The logarithmic dependence of the endpoint position on the 't Hooft cou- 
pling may signal the logarithmic branch point at infinity, which then implies 
that other complex solutions may exist, related to the original real solution 
by a non-trivial monodromy at infinity. We have studied this question nu- 
merically by solving the saddle-point equations for complex A = |A| e 1 ^ and 
slowly changing from to 2tt. We have found that the monodromy is 
actually trivial and that the saddle-point equation does not have complex 
solutions for real A. When A acquires a non-zero phase, the eigenvalue dis- 
tribution extends into the complex plane and, for small coupling, rotates in 
phase with \/A, such that for <fi going from to 2ir the eigenvalue cut makes a 
180° turn and maps to itself. The picture at strong coupling is more compli- 
cated but again the density maps to itself as soon as the phase of A changes 
by 2tt. 



4.3 Structure of the density 



When fi is large we can identify three regions in which the density has qual- 
itatively different behavior. At x ~ 1, the density is well approximated by 



the infinite-coupling solution (4.4), and varies on distances of order one. We 
call this region I. At x ~ fi, but not too close to the endpoints, the density 
changes very slowly. The scale of its variation is of order fi. There we can 
introduce the scaling variable x/fi. We call this region II. Finally, near the 
endpoints the density rapidly changes again on the scales of order one. We 
call this region III. 



to 



One can infer from (2.17) that in the region III the density is proportional 
M/A: 

A 



p(x) 



f{fi-x) 



(4.13) 



which is also consistent with the simple estimates in the previous section. 
The scaling function behaves as f(z) ~ yfz at z — > + . As we shall see later, 
f(z) ~ z~ 3 ^ 2 at z — > oo. An integral equation for f(z) of Fredholm type can 



be obtained by plugging (4.13) in (2.17) and taking /i — » oo while keeping 
H — x fixed. We will not study this equation here. 

To solve for the density in region II it is convenient to use yet another 
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integral representation for K(x): 



+00 

* , f dww coth. ttw . , 

AYx) = -/- . (4.14) 

I x — w 



The integral diverges at large w and thus this dispersion relation only holds 
after two subtractions. Namely, to make it precise we first need to shift the 
argument of K(x) by z and average over z such that the first two moments 
vanish. But this is precisely how K(x) appears in the Fredholm representa- 



tion (2.17). 



Upon substitution of (4.14) into (2.17) the double integral in y and w 



essentially acts as a unit operator as long as w 6 (— /x, //), because of the 
principal value prescription to encircle the poles. The remainder integrates 
to zero at \w\ < p,. After some transformations, we get 



■ill 



p(x) — / dy p(y) (x — y) coth ir(x — y) + x coth irx = — y p 2 — x 2 



1 f dw I fi 2 — x 2 

1 dy p(y) [(w — y) coth n (w — y) — w coth 7rw] 



(4.15) 



7r j w — x y w 2 — p? 



The advantage of this representation is that the last line is always small. For 
the major part of the eigenvalue distribution \x\ < p. Hence l/(u> — x) in 
the integrand is suppressed by a factor of order 1/p. In addition, for large x: 

£Cotli7rx ?a \x\. (4.16) 

Assuming, for definiteness, positive w and taking into account that then w—y 
is also positive, we get 



/ 



dy p(y) [(w — y) coth ir (w — y) — w coth irw] 
dyp(y) (w-y-w) = 0. 
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The last term in (4.15 ) is thus localized in a small vicinity of the distribution's 



endpoints, where p — x, w — p and p — y are all small compared to p. We can 
safely ignore this term in regions I and II, and we will also argue that this 



term is numerically small in region III. Then the equation (4.15) becomes 



p(x) — / dy p(y)(x — y) coth ir(x — y) + x coth nx 



8tt 

T 



p? — X 2 



(4.17) 



In region I, at x ~ 1, we can neglect the last term in (4.17) and solve the 



equation by the Fourier transform, which gives the asymptotic solution (4.4). 



In region II, the last term can no longer be neglected, but the density slowly 
varies with x and consequently the typical range of y's that contributes to the 
integral on the left-hand side is large, such that x — y is of order p. We can 
then use the long-range approximation for the kernel (4.16). Differentiating 
the resulting equation twice and taking into account that \x\" = 25 (x), we 
find: 



p"(x) 



- 2p(x) + 25{x) 



iixp 



A 



(p 2 - x 2 ) 



-3/2 



(4.18) 



The delta-function can be dropped, as x cannot be close to zero in region II. 
Each derivative brings in a factor of 1/p, since x scales as p, and we thus 
have: 

Anp 2 



p(x) 



A 



(p 2 -x 2 ) 



-3/2 



[X 



p). 



(4.19) 



Matching this to the solution (4.13) in region III, we find that the function 



f(z) defined there behaves at large z as 



V2n 



z -)■ oo 



(4.20) 



The solution (4.19) is non-normalizable, because it has a wrong endpoint 
behavior: (p — x)~ 3 ^ 2 instead of (p — x) 1 / 2 . To compute the normalization 
integral and thus to determine function of A we need to compute the 

density in region III and also to match the solutions in regions I and II. For 
that we need to solve eq. (4.17) for p ^> 1. For p = oo we could have used the 



Fourier transform. For large but finite interval, the solution can be obtained 
by a generalization of the Wiener-Hopf method [TBI . 
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4.4 Wiener-Hopf solution 

The Wiener-Hopf method is a generalization of the Fourier transform for the 
case when an integral equation is defined on a finite interval. The idea of 
the method is to focus on the vicinity of one end-point, and make sure that 
the boundary conditions there are correct, at the same time neglecting the 
influence of the other endpoint. This will give an accurate description of the 
density in all regions I — III for x > 0, as long as we impose the correct bound- 
ary conditions at x = //, while the solution will give a bad approximation to 
the density in the region III at x < (in the vicinity of x = —fi). But since 
the exact density is an even function of x it it sufficient to know it for positive 
x. In particular, this approximation will be sufficient for the computation of 



the Wilson loop (2.16), since it is dominated by x close to 

The Wiener-Hopf method is based on the analytic decomposition of the 
kernel in the integral equation: 

coshc 1 



2sinh 2 f G-(u)G+(u) 



where the functions 



G±(u) = . > 2 nJ 4.22 

are analytic on the upper/lower half-plane. 

The analytic properties of the Wiener-Hopf kernels are illustrated in fig. |5j 
The only singularities of G±(u) are simple poles at u — =R^ n with 

i, n = 7r(n+^J. (4.23) 



2 



The residues are 



(-2) n+1 F 2 fa + ^ 
r n = res G±(u) = 1 / ^ 2 4 (4.24) 

a/7t [n + |) T (n + 1) 



The kernels satisfy 



G±(lu) = G T (cu). (4.25) 
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o poles of C + (cv) 
• poles of C_ (cu) 

Figure 5: The analytic structure of G±(uj). 

The Wiener-Hopf solution of the integral equation, which is accurate at 
positive x, in the Fourier space is given by 



+ 



2sinh 2 ^ . 

2 F(uS) 



cosh uj cosh u 



n=0 



r„,e 



G-U)c > ~ ' (1 - F(-ix/ n )) , 

(4.26) 

where F(u) is the Fourier transform of the right-hand-side of the integral 



equation (4.15). The technical details of the derivation are collected in the 
appendix O 



The first term in (4.26) is the asymptotic solution on the infinite interval. 



The s econd term is a correction due to the linear force, which in particular 
gives (4.19) in the region II. To see this, we can replace 2sinh (u/2)/ coshcu 



by u 2 /2 — y —d 2 /2, which is justfied because the density varies very slowly in 
region II. The density then is given by the second derivative of the Wigner's 



semi-circle, in accord with (4.19). 



The last term in (4.26) subtracts the the poles of the first two terms at 
to = —%v n thus ensuring that p(u) is analytic in the lower half-plane. This 
means that the x-space density vanishes at x > fi and the solution thus 
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0.5 1.0 1.5 2.0 2.5 3.0 

Figure 6: The x-space density p(x) at A = 1000 and p = 2.4546 (solid purple 
line), compared to the infinite coupling solution Poo(x) (solid black line) and to 
Poo(x) + (p 2 — x 2 ) 3//2 , that is the region II density (dashed black line). 

satisfies the correct boundary conditions at x = fi. In fig. [6] we have plotted 
the x-spa.ce density p(x). For comparison we also plot the infinite coupling 



distribution (4.4) and the approximate solution in the region II (4.19). In 



fig. [7] the x-space density near x — // is compared to the numerical data, to 



the infinite coupling solution and to the Wigner semi-circle (4.5). 

To find fjj clS db function of A we need to impose the normalization condition 
on the density. We cannot do it by requiring that p(u = 0) = 1, because 



the solution (4.26) is approximate and gives a good approximation to the 
density only at x > 0. In particular, p(x < —p) ^ 0. p{u = 0) is given 
by the integral of the approximate solution from — oo to +oo (effectively to 
/i, since we imposed the right boundary conditions at x = p, and for x > p, 
the approximate density vanishes). Because the measure of integration is 
flat, the integral picks a contribution from the region of negative x where the 
solution does not approximate the true density with sufficient precision. 
We should instead calculate the normalization as 

fJ. +oo 

1 = 2 /^(x)= (4.27) 

J J 7T1 bj — Vo 

-oo 

using the symmetry of the exact solution under x — > —x. This trick allows 
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1.8 2.0 2.2 2.4 



Figure 7: The x-space density p(x) at A = 1000 and p = 2.4546 near x = p (solid 
purple line), compared to the infinite coupling solution poo(x) (dashed black line), 
to the Wigner semi-circle (dashed green line) and to the numerical solution of the 
saddle-point equation with N = 1200 (circles). 



us to avoid using the approximate solution beyond the region where it gives 
a controllable approximation. Plugging in the explicit expression (4.26) we 
find: 



1 + 2 



V If rf p £i(k / m~t~k / 7l) 

' m 1 n c 



n,m=0 



(F(-iv n ) - 1) 



(4.28) 



The sum over m here is dominated by the first term. The rest are exponen- 
tially small in /i. However, F(—iv n ) is exponentially large, 0( e^™), and thus 
all the terms in n should be kept. This gives the normalization condition: 



n=0 



n F( 



-IVr, 



V n + U 



2u 



(4.29) 



With the help of this condition, the density (4.26) can be rewritten in the 
form that sometimes will be more convenient: 



+ 



2 sinh 2 



cosh oj cosh u 



F(u) 



-GJu)e 



tflUJ 



oo 



r n e ^F(-iu n ) (u n - u Q ) (u - iu Q ) ^ ^ 

n=l 



{Vn + ^o) + WO) + *^n) 
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with (4.17). Then 



As a first approximation we can neglect the second line in (4.15 ) and start 

8tt 



and 



F(x) 



F(u) 



— vV - x 2 , 

87T 2 fJ,Ji(fMj) 

Xcj 



(4.31) 



(4.32) 



where J\{x) is the Bessel function of the first kind. In the lower half-plane 
we get: 

%-K 2 Hli(nv n ) e> 



F[—iv n 



XUr, 



A 



(4.33) 



Substituting this into eq. (4.29) we find that the normalization condition 



results in the equation (]4.7|) with the constant C given by 

15.16. 



ro ^ (2n + if 2 (n + 1) 



(4.34) 



The first correction due to the second line in (4.15) is computed in the 



appendix [D} Parametrically it is of the same order, but numerically turns 
out to be rather small: the correction to the normalization constant is 5C = 



—0.56. The estimate (4.8) is a combination of the leading-order result (4.34) 
and this correction. We expect that higher-order corrections are at least 
as small as the first one, which means that we know the constant C at 
a few percent accuracy. This is consistent with numerics, as discussed in 
appendix [Bj 



4.5 Wilson loop 

The Wilson loop can be computed a^] 

W(C circlc ) =p(-2m). (4.35) 



When we plug in the solution (4.26), we immediately see that the exponen- 
tially growing term in the Wilson loop vev comes from the last term in (4.26) 
which describes the solution in region III. The Wilson loop is thus determined 

4 The right-hand side is the Fourier transform of the density evaluated at uj = —2m. 
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by the behavior of the density near the endpoint of the eigenvalue distribu- 
tion, as expected on general grounds. Using (4.30) for p(uj) and (4.22) for 
G-(u), we get, omitting exponentially small terms: 

W(C circle ) -- 



-ll/r, 



2^ J>_ V- nr n e 

32^ (n + l)(3-2n) 



(4.36) 



n=l 



It is clear from (4.33) that the Wilson loop vev scales as (4.10) with 



nr r , 



^ ( n + I) 3/2 ( n + l)(3-2n) 



2.55. 



(4.37) 



The next-order correction to R is computed in appendix [D] and is not very 
big: 5R = —0.37. Altogether we get an estimate quoted in eq. (4.11). 



5 Instantons 

Since instantons can lead to a large- N phase transition [H] , we will compute 
the large- N limit of the one-instanton contribution to the partition function, 
in order to check if the moduli integration can overcome the exponential 
suppression of the instanton weight. The total one-instanton weight is given 
by 

s^2 N f n I -\2N 

^l-inst = e ^2 Re 2^ 

fc =i EI ( a k ~ a j) ( a k - a-j + 2i) 




The deformation parameters, that enter all instanton sums [13] , according to 
P] should be set to one for the M = 2 partition function on S" 4 : S\ = e 2 = 1. 
We have taken into account that the mass of the hypermultiplet is offset by 
(ei + 62) /2 = 1 [T7j, resulting in the factor of (a& + i) 2N in the numerator, 
which is the canonical contribution of 2N fundamental hypermultiplets of 
mass 1 to the one-instanton weight. 

In the large- N limit the first term in the curly brackets is exponentially 
small as long as y is not too big. Only for y of order x/iV this term becomes 
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sizable and we need to take it into account. We can thus approximate the 
integrand by expnading the exponent in the first term in 1/y: 



V 



2iV 



]1 [(y - a,) 2 + 1] 



exp 



exp 



-N 



2dj 



y y 

- (a 2 ) 



O 



y 



+ o 



The y integral becomes elementary in this approximation and we find for the 
one-instanton contribution, to the leading order in 1/JV: 



Zx- 



inst 



-4 



1 - (a?) 



Ne' 



7T 



(5.2) 



We thus conclude that the moduli integration enhances the instanton weight 
by a factor of y/N, but does not overcome the exponential suppression of the 
weight by the instanton action. 



6 Conclusions 

At weak coupling our results are in agreement with the perturbative cal- 
culation of the Wilson loop from [11]. An observation that the difference 
between the M = 4 and Af = 2 Wilson loops starts at three loops is ulti- 
mately related, as we saw, to the finiteness of the theory. The reason is the 
non-renormalization of the quadratic term in the effective action, the same 
non-renormalization property that guarantees the vanishing of the beta func- 
tion. At any order of perturbation theory, the Wilson loop vev is given by 
a combination of rational numbers and zeta-functions whose argument is 
correlated with the order of perturbation theory. It would be interesting to 
understand this transcendentality property directly from Feynman diagrams. 

At strong coupling we computed the effective string tension, which turns 
out to depend logarithmically on the 't Hooft coupling. The most straight- 
forward interpretation of this result is that the curvature of AdS§ in the dual 
geometry decreases logarithmically with A. Such an explanation implicitly 
assumes that the Wilson loop can be computed semiclassically in string the- 
ory. But the dual string theory may not have semiclassical regime^] or may be 

5 The discussion of possible quantum effects in the Wilson loop computation on the 
string side can be found in [2]. 
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semiclassical only in some approximate sense, for instance if AdS§ has small 
curvature at strong coupling, while the compact factor X 5 remains highly 
curved. If this is true, (1.3) is just an effective string tension defined through 
the parameterization (1.2) of the Wilson loop vev. In either case, it would 
be interesting to compute the Wilson loop in the dual string theory. 

Pestun's results are fairly general and potentially apply to any M = 2 
theory on S A . It would be interesting to repeat the large- N calculation of 
the Wilson loops for other superconformal, and perhaps also massive theories, 
although in the latter case the theory on S A is not equivalent to the theory 
in flat space. The closest superconformal theory to the one we studied is 
an interpolating theory with two gauge groups, which connects a Z 2 orbifold 
of M = 4 super- Yang-Mills to A^ = 2 SCYM by a continuous deformation 
[U El 12]. At the o rbifo ld point, the Wilson loop vev grows exponentially with 
\f\ according to ( 1.1 ). It would be interesting to investigate the interpolation 
to M = 2 SCYM, where the scaling is logarithmic. The transition may involve 
a non-analytic behavior [2]. 
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A Higher orders at weak coupling 

Let's approximate the function K(z) with a truncation of the Taylor expan- 



sion (3.1), i.e. 

M 



K{z) « -2 ^(-l) n C(2n + l)z 2n+1 (A.l) 



n=l 
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where M is a positive integer. Inserting the (A.l) in the expression (2.17), 
we obtain 



p(x) 



(t-;E(-i)"C(2« + D 

^ n=l 



X 



E 



2n + 1 
2A; 



n— k 



m 2 k 



E (-V2)r ^^V*" j v/^^A.2) 



fe=l v / r=0 

where we defined the r-th moment 



dz p{z)z r = (z r ) 



(A.3) 



we introduced the Pochhammer symbol (a) n = a(a — 1) . . . (a — n + 1) and 
we used 



dy 



x — 



y 



-It 



E^f (-V2H 

fc=0 



x n-l-2k ^2k 



(A.4) 



The normalization condition reads 

,2,, 2 M 



1 



/4^_ E( _ irc(2ii+i; 

^ n=l 



X 



A 

E 



2n + 1 
2A- 



^2fc At 



2+2(n-fc) 



n—k 

E 



:-i/2), 



Cn—k—r 
22{n-k-r) 



(A.5) 



where C n is the n-th Catalan number. The M moments rri2i, i = 1, . . . , M, 
can be computed using the density (A.2), giving the following system of M 
equations 



m 2i 



M 



/|_2 . .2i+2 ft 

*!li^-£(-l)»C(2» + i; 



2 2i A 



n=l 



E 

fc=i 



2n + 1 
2A; 



%—k 



m 2 kP 



2+2(i+n- 



r=0 
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where i = 1, . . . , M. The linear system (A. 6) can be be used to express the 



moments m<n in terms of /i and A and the normalization condition (A. 5) 



express \x as a function of A. The approximate solution (A. 2) therefore can 



be written explicitly at any order of approximation M and can be used to 



compute the expectation value of the circular Wilson loop (2.16). Expanding 
K{z) up to 0(z 2M+1 ), it is possible to compute /i(A) up to 0(A 3/2+M ) and the 
expectation value of the circular Wilson loop up to 0(A M+2 ). For instance, 
for M = 5 we obtain 



/X 3C(3)A 5 / 2 5C(5)A 7 / 2 7(9C(3) 2 -65C(7))A 9 / 2 
^ 2vr 256tt 5 512tt 7 65536vr 9 



+ 



3(861C(9) - 340C(3)C(5))A n / 2 



524288tt 



ii 



(-891C(3) 3 + 7900C(5) 2 + 13965C(3)C(7) - 30261C(H)) A 13 / 2 
+ 8388608tt 13 + 



(A.7) 

and 

W ^ = 1 + ^l^K9^-|^i) A3 

/ 1 27r 2 C(3)-15C(5) \ 4 

V 737280 4096tt 6 J 

f 1 3vr 4 C(3) - 65tt 2 C(5) - 12 (9C(3) 2 - 35C(7)) \ 5 

V 88473600 196608tt 8 J 
( 1 -2vr 2 C(3) + 85C(5) 

V 14863564800 7864320tt 6 
tt 2 (180C(3) 2 - 637C(7)) - 45(60C(3)C(5) - 91C(9)) \ 6 
3145728tt 10 J 
( 1 -vr 2 C(3) + 70C(5) 

V 3329438515200 377487360tt 6 

3vr 2 (108C(3) 2 - 343C(7)) - 126(110C(3)C(5) - 153C(9)) 
150994944tt 10 

27 (360C(3) 3 - 1900C(5) 2 - 3360C(3)C(7) + 4697C(11)) \ 7 
150994944tt 12 J 
+ 0(A 8 ) (A.8) 

The expression (A.7) is compared to the numerical data in in fig. [8j 
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Figure 8: The data denoted as circles are obtained by solving numerically the 
saddle-point equation with N = 100 (purple), N = 200 (green), N = 400 (red) 
and N = 600 (blue). The triangles are associated to the points extrapolated to 
N = oo. The solid line is obtained from the expression (A.7). 



B Numerical results for the endpoint posi- 
tion at strong coupling 

We have shown in section [4] that in the regime of large A, it results C = 
A e -7r ^/ 2 / y/jl where C is a constant estimated as C ~ 14.60. In order to check 
this result we have performed an extensive numerical analysis. In doing so 
we face two technical difficulties. First, it is difficult to calculate the position 
of the endpoint numerically with good precision, because the density goes 
to zero there and at a finite iV there will be few datapoints close to x — fi, 
which is visible in fig. [2j [3] and [7j We can estimate how large N should 
be. The density in the interval of interest (region III) scales as p ~ v//V^ 
(eq. ( |4.13[ )) and changes on the distances of order one. The distance between 
adjacent eigenvalues is thus of order Ax ~ l/Np ~ X/N^/JI. This should 
be at least as small as the scale of variation of the density. We thus need 
N ^ A/yZ/I ~ A(lnA) -1//2 to reach reasonable numerical accuracy. But the 
calculational cost grows as N 2 and it is impractical to go beyond N ~ 1000 
by running Mathematica on a laptop. We cannot thus reach very big A, and 
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Figure 9: The data denoted as circles are obtained by solving numerically the 
saddle-point equation with N = 400 (purple), N = 600 (blue), N = 800 (red) 
and N = 1000 (green). The triangles are associated to the points extrapolated to 
N = oo. The solid line is the interpolating function for the N = oo points and the 
dashed line is the analytical expectation C = 14.60. 

very big N, but on the other hand to estimate C we need to know /i with 
a precision that grows exponentially with fi. Our strategy was to compute 
for several values of N and then numerically extrapolate to N = oo, and 
subsequently to extrapolate in A. 

We have numerically solved the saddle-point equation in the range 50 < 
A < 5000 (as explained above we cannot get reliable results for larger A) 
considering the number of eigenvalues N = 400, N = 600, N = 800 and 
N = 1000 and computed the associated endpoints /x, see fig. [9j These data 
have been used to extrapolate the value of the endpoint fi to iV = oo for any 
given A. We have used the values of \x at N = oo to construct an interpolating 
function Cint(A) = A e ~ n ^ 2 /y/Ji such that Ci NT (oo) gives the value of the 
constant C. We obtain that C = Cint(oo) ~ 15.02, in agreement with the 
analytical result, certainly within errorbars of both calculations. 
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C Derivation of Wiener-Hopf solution 



The generalized Wiener-Hopf method replaces the Fourier transform when 
an integral equation with a difference kernel is defined on a finite interval. 
We follow the variant of the method described in the appendix of 



By Fourier transforming eq. (4.15), we get 



dui 
2^ 



cosh a; 
2 sinh 2 % 



p(u) 



1 



2 sinh' 



~F(u) 



(C.l) 



where F(u) is the Fourier transform of the right-hand side. 



If (C.l) were to hold for any x, we would immediately find but the 
equation holds only for x G (— /i, /i), and the integrand need not be equal to 
zero. The remainder can be represented as a sum of two rapidly oscillating 
functions: 



cosh to 

o — P(w 

2sinh 2 r V ; 



1 



2 sinh 



2 ui 



+ F(co) + e-^X-iu) + e^X+ico), (C.2) 



where X±(u) does not contain negative/positive frequencies and consequent- 
ly can be analytically continued to the upper/lower half-plane of complex to. 



To find p(u) we can use the decomposition (4.21) and projection opera- 



tors onto the positive and negative frequency parts. For any function, the 
projections are defined by an integral transform: 



+oo 



F{(jS) = F + (u) + J 7 - (id) 



F + (uj) 



± 



du' F{u') 



2ni uj' — u =)= iO 

— oo 

The functions F±(uj) are analytic in the upper/lower half-planes. 



(C.3) 



Multiplying both sides of (C.2) by G + (u>) e 1>1UJ , and taking the negative- 
frequency part, we find: 



p(w)e 



-1/110 



GL(w) 



G+(w) 



2 sinh' 



+ F(u) 



+ [G + (co)X4uj)e- 2i ^]_ 



(C.4) 



The functions p(u) e lflu} and F(u) e iflul are analytic in the lower half-plane, 
because their Fourier images vanish for positive x. For this reason it was 
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possible to omit the negative-frequency projection on the left-hand side. The 
function 1/ sinh 2 (a;/2) has poles in the lower half-plane, but those poles are 
cancelled by zeros of G+(u). 

Although resulting equation still contains an unknown function X_(u), 
one can argue that the last term in (C.4) is exponentially small in /i. To 
the first approximation we can neglect it. The density then is given by the 
standard solution of the Wiener-Hopf problem on a semi-infinite interval: 



1 



2 sinh' 



(C.5) 



As we have already mentioned this gives an accurate description of the den- 
sity at x > 0, and in particular satisfies the correct boundary condition at 
x — fi: p(x > fi) = 0, since p(oS) e is manifestly analytic in lower half- 
plane. But the boundary condition at x — —fi is not satisfied and we thus 
cannot trust this approximate solution at x < 0. 

The solution in the main text is obtained by closing the contour of in- 
tegration in the negative-frequency projection in the lower half-plane and 
picking up the poles of G + (u>). 



D Second iteration at strong coupling 

Here we compute the correction to the strong coupling solution due to the 



second line in (4.15): 



1 f dw / fi^ — x^ f 

SF(x) = — / 4/—= dy p(y) [(w - y) coth-K (w - y) 

n J w — x v w A — pr J 

\w\>n 

—w cotli7ru>] . (D.l) 

Assuming that x ~ fi and /i — x ~ // (this corresponds to the region II), and 
taking into account that the w integral is dominated by a small neighborhood 
of the endpoints \w ± ~ 1, we can replace l/(w — x) by l/(±/i — x). The 
whole expression then considerably simplifies: 

SF(x) « ^L= , (D.2) 

77 V z 1 — x 
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where 



.4 



dw 



w 2 — /i 2 



dy p(y) [(w — y) coth it (w — y) — w coth nw] . (D.3) 



Upon the Fourier transform, A becomes 



A 



71 



2m 2sinh 2 £ VP ; 



(D.4) 



where Hq 2 \x) is the Hankel function of the second kind. The integral can 
be calculated by closing the contour of integration into the lower half-plane 
and picking the poles at the zeros of sinh(o;/2): 



d 



71=1 



(D.5) 



k=2-ku 



Using the leading-order solution in the form (|4.30|) we find: 

A : 



d 



^ ^ di; 

n=l 

oo 



r m e ^ Um F(-iu m ) (v m - z/ ) (k + u ) 



m=l 



(u m + u Q ) (k - z/ ) (k - v m ) 



K=2irn 



At large-/i, 



/LtK 



7T 



2/iK 



and using the similar approximation (4.33) for F(—iu m ), we get: 

a 



A 



A 



(D.6) 



(D.7) 



(D.8) 



where a is a numerical constant: 

oo 

d 



mr r 



x 



1 (2m + l) 3/2 (m+1) 
r(2n+|) 



dn 2 2n n 3 / 2 r 2 (n) (2n - f) (2n - m - \) 
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1.0232. (D.9) 



This gives a correction to F(u): 



2afiJ (fj,oj) 
A 



and 



SF(— iv n 



2a/iJo(/i^r7 
A 



7TZ/„ 



fj,e 



A 



(D.10) 



(D.ll) 



Substituting this into (4.29) we find the correction to the normalization con- 



dition of the form (|4.7|) with the constant C thus shifted by 

-0.56. 



- - 2(1 y: 



^(2n + l) 1/2 (n + l) 



(D.12) 



There is no parametric suppression with respect to the leading order, but 
numerically the corrections is rather small. 

We can also calculate the correction to the coefficient R in the Wilson 



loop vev (4.10): 



nr r 



Arc 



t^i {n + l) 1/2 (n + l)(3-2n) 



-0.37. 



(D.13) 
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